Better Synchronizability Predicted by Crossed Double Cycle 
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In this brief report, we propose a network model named crossed double cycles, which are com- 
pletely symmetrical and can be considered as the extensions of nearest-neighboring lattices. The 
synchronizability, measured by eigenratio R, can be sharply enhanced by adjusting the only param- 
eter, crossed length m. The eigenratio R is shown very sensitive to the average distance L, and the 
smaller average distance will lead to better synchronizability. Furthermore, we find that, in a wide 
interval, the eigenratio R approximately obeys a power-law form as R ~ L 1 ' 5 . 
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Synchronization is observed in a variety of natural, so- 
cial, physical and biological systemsQ, and has found ap- 
plications in a variety of field including communications, 
optics, neural networks and geophysics 0, 0, El IB H 0- 
The large networks of coupled dynamical systems that 
exhibit synchronized state are subjects of great inter- 
est. In the early stage, the corresponding studies are 
restricted to either the regular networks (BQ, or the ran- 
dom ones |lOllll| . However, recent empirical studies have 
demonstrated that many real-life networks can not be 
treated as regular or random networks. The most im- 
portant two of their common statistic characteristics are 
called small-world effect^^ an d scale-free property [Tflj. 
Therefore, very recently, most of the studies about net- 
work synchronization focus on complex networks, and 
find that the networks of small-world effect and scale- 
free property may be easier to synchronize than regular 
lattices [3 El [3. 




FIG. 1: (Color online)The sketch maps of G(20,4). 
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One of the ultimate goals in studying network syn- 
chronization is to understand how the network topol- 
ogy affects the synchronizability. In the simplest case 
(see below), the network synchronizability can be well 
measured by the eigenratio R 0, 0, |3 Eof . thus 
the above question degenerates to understanding the re- 
lationship between network structure and its eigenval- 
ues. Since there are countless topological characters for 
networks, a natural question is addressed: what is the 
most important factor by which the synchroizability of 
the system is mainly determined? Some previous works 
indicated the average distance L,21\ is one of the key 
factors. However, the consistent conclusion have not 
been achieved IH |H |H El- Another extensively 
studied one is the network heterogeneity, which can be 
measured by the variance of degree distribution or be- 
tweenness distribution |26t |27| . Some detailed compar- 
isons among various networks have been done, indicating 
the network s ync hronizability will be better with smaller 
heterogeneity [3 E3 ■ However, a well-known coun- 
terexample is the regular networks with homogeneous 
structure while displaying very poor synchronizability. 
Because of the networks used for comparing in previ- 
ous studies are of both varying average distances and 
degree variances, the strict and clear conclusions can 
not be achieved. In addition, some researchers deem 
that the more intrinsic ingredient leading to better syn- 
chronizability is the randomicitv|30j . that is to say, the 
intrinsic reason making small-world and scale-free net- 
works having better synchronizability than regular ones 
is their random structures. Therefore, if one wants to 
clearly show how L affects the network synchronizability, 
he should investigate the networks of different L but the 
same degree variance. And if he wants to assert that it is 
not the randomicity but smaller (or longer) L resulting in 
the better synchronizability, the deterministic networks 
are required. 

In this brief report, we proposed a deterministic net- 
work model named Crossed Double Cycles (CDCs for 
short). The CDCs are of degree variance equal to zero, 
and by adjusting the only parameter m, named crossed 
length, the average distance of CDCs can be changed. By 
using this ideal model, we demonstrate that the smaller L 
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FIG. 2: (Color online) The average distance of the CDCs. The 
black squares, red circles, blue triangles and green pentagons 
represent the case of N =1000, 2000, 3000 and 4000, respec- 
tively. The smaller solid symbols and larger hollow symbols 
represent the simulation and analytic results, respectively. 



FIG. 3: (Color online) Eigenvalues. The network size N = 
4000 if fixed, and the black squares, red circles and green tri- 
angles denote the numerical results of eigenvalues for the cases 
m = 2, m = 3, and m = 4, respectively. The corresponding 
curves represent the analytical solutions. 



will result in better synchronizability, and provide a use- 
ful method to enhance the synchronizability of nearest- 
neighbor coupling networks. 

In network language[3l|, the cycle Cm denotes a net- 
work consisting of N vertices x\,x%,- ■■ , xjy. These N 
vertices is arranged as a ring, and the nearest two vertices 
are connected. Hence, Cm has ./V edges connecting the 
vertices Xix 2l x 2 x%, • • • , xm-iXn, and Xn%i- The CDCs, 
denoted by G(N,m), can be constructed by adding two 
edges, called crossed edges, to each vertex in Cat. The 
two vertices connecting by a crossed edge are of distance 
to in Cat. For example, the network G(N, 3) can be con- 
structed from Cm by connecting X1X4, X2X5, • • • , xm-±X2, 
and x mx$ ■ And the network G(N, 2) is isomorphic 32] to 
a one-dimensional lattice with periodic boundary con- 
ditions wherein each vertex connects to its nearest and 
next-nearest neighbors. A sketch map of G(20, 4) is 
shown in Fig. 1. 

Clearly, all the vertices in G(N, m) are of degree 4, thus 
the degree variance is equal to 0. Furthermore, G(N, m) 
is vertex-transitivity, that is to say, for any two vertices 
x and y in G(N, m), there exists an automorphism map- 
ping 6 : V(G) -> V(G) such that y = 6{x). The vertex- 
transitivity networks are completely symmetric, which 
are of particular practicability in the design of topolog- 
ical structures of data memory allocation and multiple 
processor systems [33) . 

Denote L(k) the average distance of Ck+i, we have 



L(k) 




k is even 
k is odd 



(1) 
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For N 3> to, we assume N can be exactly divided by to 
and denote k = — . Since G(N,m) is vertex-transitivity, 
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FIG. 4: (Color online) The smallest nonzero eigenvalue 71 (#1) 
vs m. The black squares, red circles and green triangles de- 
note the cases N = 10000, N = 250000, and N = 1000000, 
respectively. The solid lines are of slope 2 for comparison. 
The inset shows the cut-off point m c as a function of network 
size. 



the average distance of G(N, m) is equal to the aver- 
age distance between vertex x\ to all other vertices. 
The network G(N, m) contains k end-to-end C m +i as 

X\X2 ■ ■ ■ X m+ i,X m+ iX m+ 2 ■ ■ ■ X 2m+ i, • • • , XN+l-m ' ' ' XffXl. 

Then, go from the vertex x\ to a certain vertex Xi can 
be divided to two processes. Firstly, travel through the 
crossed edges to the nearest vertex that belongs to ar^'s 
cycle mentioned above. Secondly, pass by a shortest 
path restricted in this cycle to X{. For example, the 
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FIG. 5: (Color online) R vs L. The black squares, red cir- 
cles, blue triangles and green pentagons represent the case of 
N =1000, 2000, 3000 and 4000, respectively. The inset shows 
the same data in log- log plot, indicating that the eigenratio 
R approximately obeys a power-law form as R ~ L 1 ' 5 . The 
solid line is of slope 1.5 for comparison. 



path from x\ to xio in G(20,4) is x\ — ► x§ — > xg — > ojio. 
The first two edges are the crossed edges and identified 
by red lines in Fig. 1, and the last edge is in the cycle 
XgXioXiia^ccis. Hence one can obtain the average 
distance of G(N, m) for N ^S> m as: 



L G (N, m) = L{m) + L(k) - 1. 



(2) 



Figure 2 shows the simulation results of Lo(N,m) for 
N = 1000,2000,3000,4000 and m < 100, which accu- 
rately agree with the analytic ones. 

In succession, we investigate the changes of CDCs' 
synchronizability with m. Consider N identical dynam- 
ical systems (oscillators) with the same output func- 
tion, which are located on the vertices of a network and 
coupled linearly and symmetrically with neighbors con- 
nected by edges of the network. The coupling fashion 
ensures the synchronization manifold an invariant man- 
ifold, and the dynamics can be locally linearized near 
the synchronous state. The state of the ith oscillator 
is described by x% and the set of equations of motion 
governing the dynamics of the N coupled oscillators is 



JY 



(3) 



where x 1 = F(x l ) governs the dynamics of individual 
oscillator, H(x J ) is the output function and a is the cou- 
pling strength. The N x N Laplacian G is given by 
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ki for i = j 
- 1 for j e A., . 
otherwise 



(4) 



Because of the positive semidefinite of G, all the eigen- 
values of it are nonnegative reals and the smallest eigen- 
value #o is always zero, for the rows of G have zero sum. 
Thus, the eigenvalues can be ranked as 9q < 9\ < ■ ■ ■ < 
On— i- The ratio of the maximum eigenvalue 9n-i to the 
smallest nonzero one 9% is widely used to measure the 
synchronizability of the net work |!7lll9j|. if the eigenratio 
R = 0n-i/9i satisfies 



R < 0.2/0.1, 



(5) 



we say the network is synchronizable. The right-hand 
side of this inequality depends only on the dynamics of 
individual oscillator and the output function|l9j. while 
the eigenratio R depends only on the Laplacian G. R in- 
dicates the synchronizability of the network, the smaller 
it is the better synchronizability and vice versa. In this 
brief report, for universality, we will not address a par- 
ticular dynamical system, but concentrate on how the 
network topology affects eigenratio R. 

Since the Laplacian for any CDC is shift invariant, 
the eigenvalues can be calculated from a discrete Fourier 
transform of a row of the Laplacian matrix |20l| . Denote 
7i (i = 0, 1, • • • , N — 1) the ith eigenvalue |34j, it reads 



2iri 2nim 

7, = 2 2 — cos cos 

1 N N 



(6) 



Figure 3 shows the numerical results of eigenvalues, 
which accurately agree with the analytical solutions. 
Clearly, for odd m and even N, the maximal eigenvalue is 
9n—i — In/2 = 8. However, if m is even or N is odd, the 
maximal eigenvalue 9n-i is smaller than 8 and can not 
be expressed sententiously. The smallest nonzero eigen- 
value 9\ equals 71, and can be approximately obtained 
under the condition N 3> m 



7i 



— (1 



m 2 ). 



(7) 



Figure 4 reports 71 as a function of m, which scales as m 2 , 
as predicted by Eq. (7), before reaching a cut-off point 
m c . The numerical value of m c is also shown in the inset 
of Fig. 4, which accurately obeys the form m c — \^N. 

Because of the cut-off in 9\, the fluctuations in 9m-i for 
even m, and the relatively complex relationship between 
m and L, we can not obtain a straightforward expres- 
sion to comprehensively depict the relationship between 
R and L. In Fig. 5, we only report the numerical results 
about how the average distance affects the network syn- 
chronizability. One can see clearly, the network synchro- 
nizability is very sensitive to the average distance; as the 
increase of L, the eigenratio R sharply spans more than 
three magnitudes. And the network synchronizability is 
remarkably enhanced by reducing L. When the crossed 
length m is not very small or very large (comparing with 
N), the networks with the same average distance have ap- 
proximately the same synchronizability no matter what 
the network size is. More interesting, the calculated re- 
sults indicate that the eigenratio R approximately obeys 



4 



a power-law form as R ~ i 1,5 in a wide interval of L (see 
the inset of Fig. 5). 

To sum up, we propose an ideal network model, and in- 
vestigate its synchronizability. The results indicate that 
the average distance is an important factor affecting the 
network synchronizability greatly. The smaller average 
distance will lead to better synchronizability. This is 
similar to the communication systems, wherein the aver- 
age distance is one of the most important parameters to 
measure the transmission delay (or time delay) encoun- 
tered by a message travelling through the network from 
its source to destination, and the smaller average distance 
means higher efficiency for homogeneous networks. Very 
recently, by numerical studies, some authors think that 
there may exist some common features between dynam- 
ics on communication networks (tra ffic and diffusion) and 
network synchronization }^ I35T.I361 137| . Since in the for- 
mer dynamics shorter L will lead to greater throughput 



and fast spread, the underlying common features pro- 
vide a possible explanation why shorter average distance 
corresponding to better synchronizability. 

The CDCs are natural extensions of the lattice of near- 
est neighbors, they are symmetric and with better syn- 
chronizability, thus have great potential in the applica- 
tions for designing of topological structures of distributed 
processing systems, local area networks, data memory al- 
location and data alignment in single instruction multi- 
ple data processors 33J. In fact, the processor network 
of one kind of the earliest parallel processing comput- 
ers is G(16, 4)|33|. Besides in synchronization, the cross 
method also has been applied in communication sys- 
tems. For example, the crossed cubes have much larger 
throughput than hypercubes thus are widely used in de- 
signing parallel computing networks 39, 40]. 
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